

Td = dindgen(1001)*0.03d + 10.d

print,m4_stat(Td)

gamma = FLUX2TAU_BGPS( Td )


myps,'./irdc_dist_model/analysis_plots/plot_gamma_bgps.eps',xsize=5

;multiplot,[2,1],xgap=0.04

plot,Td,gamma,xtit='Dust Temperature [K]',$
     ytit=cgSymbol('alpha')+'(T!dd!n) [Jy!u-1!n]',thick=5


;; multiplot,/doyaxis


;; plot,[0],[0],xtit='Dust Temperature [K]',ytit='f!ddata!n',thick=5,/nodata,$
;;      xr=[10,40],yr=[-1,1]

;; vline,/horiz,0,linestyle=1


;; C = [0.1,0.25,0.4]
;; S = [0.1,0.3,1.0]

;; colors = cgColor(['BLU5','GRN4','RED5'])
;; linest = [0,2,3]

;; FOR j=0,2 DO BEGIN
;;    FOR k=0,2 DO BEGIN
      
;;       fdata = ( 1.d - ( C[j] / (1.d - exp(-S[k] * gamma) ) ) ) ;> 0.d
      
;; ;; plot,Td,fdata,xtit='Dust Temperature [K]',ytit='f!ddata!n',thick=5
      
;;       oplot,Td,fdata,thick=5,linestyle=linest[j],color=colors[k]
;;    ENDFOR
;; ENDFOR



;; multiplot,/reset
;; multiplot,/default
;; multiplot,/reset

myps,/done

END
